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FOREWORD 


This  report  describes  a  portion  of  the  results  obtained  on  NASA  Grant 
NSG  3044.  This  work  was  done  under  subcontract  to  the  University  of  Illinois, 
Urbana,  with  Prof.  S.S.  Wang  as  the  Principal  Investigator.  The  prime 
grantee  was  the  Massachusetts  Institute  of  Technology,  with  Prof.  F.J.  McGarry 
as  the  Principal  Investigator  and  Dr.  J.F.  Mandell  as  a  major  participant. 

The  NASA  -  LeRC  Project  Manager  was  Dr.  C.C.  Chamis. 

Efforts  in  this  project  are  primarily  directed  towards  the  development 
for  finite  element  analyses  for  the  study  of  flaw  growth  and  fracture  of 
fiber  composites.  This  report  presents  study  of  boundary  layer  problems  in 
composites  based  on  an  exact  solution. 
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ABSTRACT 


A  study  of  boundary-layer  stress  singularities  in  multilayered  fiber- 
reinforced  composite  laminates  is  presented.  Based  on  Lekhnitskii ' s  stress 
potentials  and  the  theory  of  anisotropic  elasticity,  formulation  of  the  problem 
leads  to  a  pair  of  coupled  governing  partial  differential  equations.  An  eigen¬ 
function  expansion  method  is  developed  to  obtain  the  homogeneous  solution  for 
the  governing  P.D.E's.  The  order  or  strength  of  boundary-layer  stress  singu¬ 
larities  is  determined  by  solving  the  transcendental  characteristic  equation 
obtained  from  the  homogeneous  solution  for  the  problem.  Numerical  examples  of 
the  singular  strength  (or  singular  eigenvalues)  of  boundary-layer  stresses  are 
given  for  angle— ply  and  cross— ply  composites  as  well  as  the  cases  of  more  general 
composite  lamination. 


ii 


1 .  INTRODUCTION 


The  response  of  a  multilayered  fiber-reinforced  composite  laminate  near 
its  geometric  boundaries  has  been  a  subject  of  intensive  investigation  during 
the  last  decade.  Both  experimental  studies  and  approximate  analytical  solu¬ 
tions  have  indicated  that  complex  stress  states  with  rapid  change  of 
gradients  occur  along  the  edges  of  composite  laminates,  for  example,  Refs  [1- 
15].  This  phenomenon  is  considered  to  result  from  the  presence  and  inter¬ 
actions  of  geometric  discontinuities  of  the  composite  and  materials  discon¬ 
tinuities  through  the  laminate  thickness .  The  anomaly  has  been  found  to  occur 
only  within  very  local  region  near  the  geometric  boundaries  of  a  composite 
laminate,  and  is,  therefore,  frequently  referred  to  as  "boundary-layer  effect" 
or  "free-edge  effect"  —  a  problem  unique  to  composite  laminates  and  not 
observed  in  homogeneous  solids  in  general.  It  has  been  shown  further  that 
the  boundary-layer  effect  is  three-dimensional  in  nature  and  not  predictable 
by  the  classical  laminate  theory  (CLT)  [16,17].  The  boundary-layer  effect  is 
apparently  one  of  the  most  fundamental  and  important  problems  in  the  mechanics 
and  mechanical  behavior  of  composite  laminates.  The  high  stresses  developed 
in  the  boundary- layer  region  coupled  with  the  low  interlaminar  strength  are 
certainly  of  critical  significance  in  aggravating  the  failure  of  composite 
materials  and  structures.  For  example,  boundary-layer  stresses  have  been 
observed  to  be  responsible  for  the  initiation  and  growth  of  local  hetero¬ 
geneous  damage  in  the  forms  of  interlaminar  (delamination)  and  intralaminar 
Ctransverse  cracking)  fracture  in  composite  laminates  under  static  loading 
[3,18].  They  are  considered  to  have  even  greater  effects  on  the  long  term 
strength  of  composite  laminates  under  cyclic  fatigue  loading  [19,20]. 
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While  the  significance  of  boundary-layer  effects  has  long  been  recognized, 
research  progress  on  this  subject  has  been  relatively  slow.  The  situation  is 
apparently  caused  by  the  inherent  complexities  involved  in  the  problem:  the 
strong  anisotropy  of  mechanical  properties  of  each  individual  ply,  the  abrupt 
change  of  materials  properties  through  the  laminate  thickness ,  the  geometric 
discontinuity  along  laminate  boundaries,  and  the  coupling  between  in-plane 
and  transverse  deformations  and  stresses  near  the  edges  of  the  composite  lami¬ 
nate.  According  to  Pagano  [14],  analytical  studies  to  date  may  be  roughly 
classified  into  two  general  categories:  approximate  theories  and  numerical 
solutions.  The  first  approximate  solution  for  finite-width  composite  laminates 
was  proposed  by  Puppo,  et  al.  [4]  based  on  a  laminate  model  containing  aniso¬ 
tropic  laminae  and  isotropic  shear  layers  with  the  interlaminar  normal  stress 
being  neglected  throughout  the  laminate.  Other  approximate  theories  were  also 
attempted  to  examine  the  problem  such  as  the  extension  of  the  higher-order 
plate  theory  [21]  by  Pagano  [10],  the  perturbation  method  by  Hsu,  et  al.  [12], 
and  a  boundary-layer  theory  by  Tang,  et  al.  [11].  Recently,  Pagano  [14,15] 
has  developed  an  approximate  theory  based  on  assumed  in-plane  stresses  and 
the  use  of  Reissner’s  variational  principle.  Even  though  there  is  no  singu¬ 
larity  involved  in  the  formulation,  the  approach  has  certain  features  signifi¬ 
cantly  important  in  objectively  determining  detailed  laminate  stress  fields. 

The  study  of  edge  stresses  in  composites  by  using  a  numerical (finite  differ¬ 
ence  method  was  apparently  first  made  by  Pipes,  et  al.[5].  Isakson  and 
Levy  [6]  developed  a  finite-element  scheme  containing  membrane  elements, 
which  closely  resemble  the  laminate  model  of  Puppo  et  al.  [4].  Later  finite 
element  studies  on  this  subject  by  Wang,  et  al.  [13]  and  Herakovich,  et  al. 

[18]  led  to  numerical  solutions  similar  to  that  given  by  Pipes  [5].  Due  to 
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the  singular  nature  of  the  problem,  a  large  number  of  elements,  especially 
through  the  thickness  direction,  are  required  in  conjunction  with  a  lengthy 
extrapolation  procedure  in  order  to  achieve  satisfactory  solutions  even  for 
a  simple  two  or  three  layer  laminate.  Improved  finite-element  methods  by 
using  a  more  complex  element  stiffness  formulation  based  upon  Maxwell  stress 
functions  [7]  and  by  hybrid-stress  elements  [22]  have  been  able  to  achieve  an 
expedient  computation  with  significantly  less  elements.  Unfortunately,  the 
refinements  do  not  guarantee  [23]  the  convergence  and  accuracy  of  the  numer¬ 
ical  solutions  because  of  the  singular  nature  of  the  boundary-layer  stress 
field.  That  is,  with  each  more  refined  analysis,  numerical  values  of  the 
maximum  interlaminar  stresses  are  shown  to  rise  with  continuously  decreasing 
element  size.  The  quest,  apparently,  is  to  show  that  a  stress  singularity 
exists  at  the  edge  of  a  composite  laminate. 

From  a  linear  elasticity  point  of  view,  it  is  well  known  that  stress 
singularities  are  prevalent  at  the  comers  of  geometric  boundaries  joining 
dissimilar  materials  (see,  for  example,  [24-26]).  Unfortunately,  the  search 
for  the  order  of  stress  singularities  for  the  boundary-layer  region  in  a  com¬ 
posite  laminate  containing  anisotropic  plies  has  not  been  successful  to  date, 
to  the  authors'  knowledge.  Since  the  singular  boundary-layer  stresses  are 
extremely  localized  in  nature,  the  precise  nature  of  the  boundary-layer  effect 
will  not  be  fully  understood  until  the  exact  order  of  the  stress  singularities 
is  defined.  In  this  paper,  the  first  in  succession,  a  rigorous  theoretical 
investigation  of  the  free-edge  stress  singularity  in  composite  laminates  is 
presented. 

In  the  next  section,  a  mathematical  model  and  basic  equations  for  each 
lamina  of  the  composite  are  presented.  Based  on  the  theory  of  anisotropic 
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elasticity  and  Lekhnitskii's  stress  potentials  [27],  a  pair  of  linear 
governing  partial  differential  equations  is  derived.  Appropriate 

near  field  boundary  conditions,  end  loading  conditions  and  interface  contin¬ 
uity  conditions  are  given  also.  The  homogeneous  solution  for  the  problem  is 
obtained  in  Section  3  by  an  eigenfunction  expansion  method.  The  solution  proce¬ 
dure  used  to  evaluate  the  exact  order  of  the  boundary-layer  stress  singularity 
is  presented.  Degenerated  cases  of  commonly  used  cross-ply  composite  laminates 
are  examined  also.  Numerical  examples  of  determining  the  edge  stress  singu¬ 
larities  for  graphite-epoxy  composite  laminates  with  various  fiber  orientations 
are  given  in  Section  4.  As  will  be  shown  later,  the  existence  of  the  free- 
edge  stress  singularity  in  composite  laminates  is  proven  mathematically  in 
this  paper.  It  settles,  once  and  for  all,  the  previous  conjecture  of  boundary- 
layer  stress  singularities  in  composite  materials,  and  provides  a  rigorous 
mathematical  method  for  determining  the  exact  value  of  the  edge  stress  singu¬ 
larity,  which  is  the  fundamental  basis  for  the  boundary  layer  theory  in  com¬ 
posite  materials  and  structures. 
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2 .  FORMULATION 


2.1  Basic  Equations 

Consider  a  composite  laminate  composed  of  fiber-reinforced  plies  with 
constitutive  equations  described  by  generalized  Hooke’s  law  in  the  con¬ 
tracted  notation  as 


e. 

1 


S.  .0. 

13  3 


(i,j  =  1,2, 3, 4, 5, 6),  (1) 


where  the  repeated  subscript  indicates  summation  and  is  the  compliance 
tensor.  The  engineering  strains,  e^,  in  Eq  1  are  defined  in  a  Cartesian  co¬ 
ordinate  system  by 


e,  = 


e,  = 


e  = 

X 


yz 


9u 

9v 

9w 

3x’ 

2  y  9y 

^3  ■  ^z  "  9z’ 

9w  L 

9v 

9w  , 

9u 

9u 

'  3y 

9z’ 

"5  ■  ''xz '  sir  + 

9z’ 

=  Y  - 

6  xy  oy 

+  & 


where  u,  v  and  w  are  the  components  of  displacements.  The  stresses,  o^,  are 
defined  in  an  analogous  manner  in  the  Cartesian  coordinate  system. 

The  composite  laminate  considered  here  has  a  finite  width  and  is  sub¬ 
jected  to  surface  tractions  acting  in  planes  normal  to  the  generator  of 
the  lateral  surface  and  not  varying  along  the  generator,  i.e.,  along  the  z- 
axis  (Fig.  1).  The  composite  is  assumed  to  be  sufficiently  long  that,  in 
the  region  far  from  the  ends,  the  end  effect  is  neglected  by  virtue  of  the 
Saint  Venant  principle.  Consequently,  the  stresses  in  the  laminate  are  in¬ 
dependent  of  the  z-coordinate .  The  case  of  a  finite-width  composite  laminate 
subjected  to  a  uniform  axial  strain,  z-axis  has  been  inten¬ 

sively  studied  by  many  researchers  [5-13  ]  .  The  special  case  in  which 
stresses  and  displacements  are  independent  of  z  corresponds  to  the  well  known 
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generalized  plane  deformation  [27] .  Under  these  assumptions,  the  equations  of 
equilibrium  without  body  force  read 


3a  3t 

+  . xy,  ^ 

3x  3y 


3t  3a 

3x  3y 


Following  the  procedure  in  [27],  it  can  be  shown  after  some  mathematical 
manipulation  that  the  general  expressions  for  displacements  and  the  stress 
component  a^  have  the  following  form: 

1  5 

u  =  -  y  “  A^j^yz  +  U(x,y)  +  -  u^y  +  u^,  (4a) 

1  2 

V  =  -  —  A^xz  +  V(x,y)  +  w^x  -  w^z  +  v^,  (4b) 

w  =  (A  x  +  A„y  +  A-)S»-Z  +  W(x,y)  +  U3^y  -  oj„x  +  w  ,  (4c) 

X  Z  J  JO  i  /  o 


",  -  V"  V  b  - 


(j  =  1,2, 4, 5, 6) 


The  unknown  functions,  U,  V  and  W,  depend  on  x  and  y  only,  and  can  be  shown 
easily  to  obey  the  following  relations: 

3TT 

11=32^0^  +S23(AiX  +  A2y  +  A3),  (5b) 

^53  +  V- 

I7  ■  'hfl  *  *2?  *3'  '  <=<*> 

3TI  9V  ~ 

37  ^63  ^^1^  ^2^ 


where 


S..  =  S..  -  S.„S.„/S„„, 
13  ij  i3  33  33’ 


(i,3  =  1,2,4,5,6) . 


It  is  obvious  that  the  constants,  u  ,  v  ,  w  and  u .  (i  =  1,2,3)  in  Eqs  4a-4d 

O  O  O  1  77/  -I 

characterize  the  rigid-body  translations  and  rotations  of  the  solid.  A^  and 
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represent  the  bending  of  the  laminate  in  the  x-z  and  y-z  planes. 
characterizes  the  uniform  axial  extension  of  the  composite  laminate,  and 
A^,  the  relative  angle  of  rotation  about  the  z-axis. 

2.2  Governing  Partial  Differential  Equations 

Introducing  Lekhnitskii's  stress  potentials,  F(x,y)  and  'i'(x,y)  [27], 
such  that 


3x3y’ 


one  can  show  that  the  equations  of  equilibrium  are  satisfied  identically. 
Eliminating  U  and  V  from  Eqs  5a,  5b  and  5e  and  W  from  Eqs  5c  and  5d  by  dif¬ 
ferentiation,  we  obtain  t^ie  following  system  of  governing  partial  differential 
equations  for  the  problem: 


[  L3F  +  L2'i'  -  -  2A^  +  A^S^^  -  A2S33, 
< 

L.F  +  L_¥  =  0, 

4  j 


where  and  are  linear  differential  operators  defined  as 

^2  "  ^44  "^^45  ax8y  ^55  iy^  ’ 

^3  "  ■  ^24  ^  ^^25  ^46^3!^  “  ^^14  ^56^3x3y^  ^15  3y^’ 


\  "  ^22  ^  ■  ^^26  sAy  ^^^12  ^66^3x^3y^  ^^16  3x3 y^  ^11  ^ 


8 


2.3  Boundary  and  End  Conditions 

Assuming  that  the  edges  of  a  composite  laminate,  are  traction  free 

and  the  interface  of  the  mth  and  (OT+l)th  plies  is  a  straight  line  meeting 
the  traction-free  edge  at  a  right  angle  (Fig.  1),  one  can  obtain  the  following 
boundary  conditions  along  36^: 

o  =  T  =  T  =  0.  (8) 

X  xy  X2 

The  conditions  at  the  ends  of  the  composite  laminate  may  have  the  form  from 

the  statically  equivalent  loads  as 

’  r  f  c  r  r 

T  dx  dy  =  0,  T  dx  dy  =  0,  dx  dy  =  P^, 

JJgXz  jj^yz  jjg 

*  r  r  r  r  r 

a  y  dx  dy  =  M  ,  a  x  dx  dy  =  M  ,  ("t,,  x  -  t  y)dx  dy  =  M  , 

jJgZ  xJJgZ  yjjgyz  ^ 

(9) 

where  the  integrals  are  taken  over  the  entire  area  B  of  the  cross  section, 

and  P  ,  M  ,  M  and  M  are  the  applied  force,  bending  moments  and  twisting 
z  X  y  t 

moment  acting  on  the  ends,  respectively. 

2.4  Interface  Continuity  Conditions 

Consider  a  portion  of  the  laminate  cross  section  composed  of  the  mth  and 
(m+l)th  fiber-reinforced  laminae,  as  shown  in  Fig.  1.  Assuming  that  the 
plies  are  perfectly  bonded  along  the  interface  SB^.,  one  can  immediately  es¬ 
tablish  the  continuity  conditions  of  the  stresses  and  displacements  along 
the  interface  as  the  following: 


,(m) 

n(m) 

+ 

(m+1) 

*  -a 

_  ^Cm+l) 

(m+1) 

(lOa) 

X 

X 

xy 

y 

X 

X 

xy 

y 

(m) 

(m) 

XI 

I.W 

(m+1) 

n 

- 

(m+l) 

^ .  > 

(10b) 

xy 

X 

y 

y 

xy 

X 

y 

y 

(10c) 


and 


(m)  (m)  ,  (m)  (jn)  (m+l)  (m+1)  (w+1)  (m+1) 

T^  n  +T  n  =-T  n  -t  n 

xz  X  yz  y  xz  x  yz  y 

(m)  (m+1)  (m)  (m+1)  (m)  _  (m+1) 

u=u  ,  v=v  ,  w-w  , 


(lOd-f) 


where  the  superscripts  denote  the  mth  and  (m+l)th  plies  in  a  composite 

laminate,  and  n  and  n  are  components  of  unit  outward  normal  to  the  inter- 
X  y 

face. 
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3.  HOMOGENEOUS  SOLUTION  AND  FREE-EDGE  STRESS  SINGULARITY 

The  governing  equations,  7a  and  7b,  are  coupled,  linear  partial  differen¬ 
tial  equations  with  constant  coefficients  related  to  the  anisotropic  elastic 
constants  of  each  individual  lamina.  With  the  aid  of  aforementioned  near-field 
boundary  conditions  and  interface  continuity  conditions,  the  homogeneous  solu¬ 
tion  for  the  governing  P.D.E.'s  can  be  determined  easily.  The  homogeneous 
boundary  conditions  and  interface  continuity  conditions  also  provide  the  in¬ 
formation  for  determining  the  important  strength  or  order  of  the  free-edge 
stress  singularity  in  a  composite  laminate,  which  is  the  major  concern  in 
this  paper. 

According  to  Lekhnitskii  [27],  the  homogeneous  solution  for  the  governing 
partial  differential  equations  has  the  general  form  as 


F(x,y)  =  I  F  (x  +  p  y). 


’P(x,y)  =  I  (lla-b) 


=1  k=l 

where  the  prime  (')  in  Eq  11b  denotes  differentiation  of  the  function 

F^(x  -I-  with  respect  to  its  argument,  and  the  coefficients  are  the  roots 

of  the  following  algebraic  characteristic  equation 


where 


Cy)  - 

(y)  =  0, 

(12a) 

\  =  - 

2^^k^ 

-  '‘4(Uk)/»3(Pk>' 

(12b) 

S55P^  - 

2S^5y  + 

(12c) 

^-3(y)  = 

hs’'"  - 

^46^^  ^24’ 

(12d) 

^^Cy)  = 

hi'*-  - 

+  (25^2  +  i66)p' 

-  2S2gy  +  5^2 • 

(12e) 
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It  can  be  shown  that  Eq  12a  cannot  have  a  real  root  (thus,  have  to  appear  as 
complex  conjugates) ,  and  are  analytic  functions  of  the  complex  variables 

=  X  +  y^y  =  r(e^^  +  X^e  ^^)/(l  +  X^)  with  X^  =  (1  +  (1  -  and 

r  and  9  being  components , of  polar  coordinates.  Substituting  the  expressions 
of  F(x,y)  and  'i'(x,y),  Eqs  11a  and  11b,  into  Eqs  6a-6e,  the  homogeneous  solu¬ 
tions  for  stresses  may  be  expressed  in  terms  of 


4'”  =  I 


(13a) 


(13b) 


-  -I,  \ 


(13c) 


“k\ 

k=l 


(13d) 


(13e) 


The  expressions  for  displacement  components  may  be  obtained  directly  from 
Eqs  5,  7  and  13  with  omission  of  the  terms  which  are  to  be  included  in  the 
particular  solution. 


u'**’ -  I  P.  <(V' 


k  k'  k" 


(14a) 


-  I 


(14b) 


.(i^)  = 


(14c) 
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where 


’k 

'k 


ll^k  + 

^12  ~  ^14'^k  ■'*  ^15^k^k  ■  ^b’^k’ 

(14d) 

IZ^k  ■" 

^22^^k  ■  ^24^k^^k  ■*■  ^25’^k  ■  ^26’ 

(14e) 

U^k^ 

^24^’^k  “  ^44V’^k  ■*■  ^45\  "  ^46' 

(14f) 

We  now  choose  the  form  of  as 


,6+2 


\(Zk)  =  Vk 


(15) 


where  and  6  are  arbitrary  complex  constants  to  be  determined  later.  Sub¬ 
stituting  Eq  15  into  Eqs  13  and  14  gives 


“f’  =  I  uC  rj . 

k=l 


2„5 


-2  7$. 

■'k^k“k  ■  '^k+S  ^k 


(16a) 


=  I  1=^4  +  c.  Z?1 , 


"yS'  ■  :l  *  Vs  ’ 

•'  k=l 


k=l 

3 


Z^' 

"k“k  '  '^k+3  k' 


(16b) 


(16c) 


T =  Z  [C^n^p^z^  +  p^  Z^I . 

k=l 


—  —  ■ 

k"k^k“k  ■  "k+3  "^k  ’^k 


(16d) 


,(h)  _  ?  „5 


xy 


k=l 


=  -I  tVk^k  '=k«  “k  ’ 


(16e) 


and 


■  I  'Wk"  +  <=k*s  Pv  +  1). 


.ft)  _ 


k=l 

3 

I 

k=l 

3 


k+3  ‘^k  k 


I  'wr^  ■"  '>k  + «• 


k+3  ^k  k 


k=l 


(17a) 


(17b) 


(17c) 
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where  the  overbar  denotes  the  complex  conjugate  of  the  associate  variable. 

For  convenience,  we  drop  the  superscript  h  associated  with  the  above  homogeneous 
solutions  for  stresses  and  displacements  in  this  paper. 


The  homogeneous  solutions  are  required  to  satisfy  the  homogeneous  boundary 
conditions  and  interface  continuity  conditions.  This  leads  to  a  standard 
eigenvalue  problem  for  determining  the  values  of  d .  It  is  noted  that  5  general¬ 
ly  appears  as  a  set  of  complex  conjugates,  which  enable  to  make  Eqs  16  and  17 
real  functions  by  superposition.  Furthermore,  the  value  of  5  is  required  to 
satisfy  the  condition 

Re[5]  >  -1  (18) 

to  ensure  the  finiteness  of  displacement  components  at  the  origin,  where  Re 
represents  the  real  part  of  5 . 

To  expedite  further  developments,  we  transform  the  stress  and  displace¬ 
ment  components  from  Cartesian  coordinates  to  polar  coordinates.  Thus,  we 


have 


"ee  =  J,  ^  S+3  "ik  4>’ 

k=l 


(19a) 


^ez  "  J.  ^^k®2k\  ^k+3  ®2k  ’ 

k=l 


(19b) 


"^er  "  ^^k®3k\  ^k+3  ®3k  V’ 

k=l 


(19c) 


a  =  y  (c,  H,.  Z,  +  C-  H,,  Z  ), 
rr  ,  ^ ,  k  4k  k  k+3  4k  k 
k=l 


(19d) 


^rz  “  ^V5k\  ^k-l-3  ^Sk  ’ 

k=l 


(19e) 
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and 


15 


^  ^k  '  2^  ^“k  '  ^k+3  jk  ^  2'  ^"k  ^  2^ 

e=l  -' 


(22b) 


I  +  c«  r^f  1  -  -  o,  (22c) 

k  rk  k+3  rk  k  rk  k+3  rk 


k  rk 


k+3  rk 


(j  =  1,2,3;  r  =  1,2,3,4,5,6) , 


where  H  ,  (t)  and  H.,  (— ^)  are  values  of  H. .  evaluated  at  0  =  ^  and  9  -  -  -r, 
jk  2  jk  2'  ij  2  / 

respectively;  f2^(9)  are  defined  as 


f2j^(9)  =  (e^®  +  l^e~^®)/(l  +  X^)  , 


^Ik  “  ^2k  ”  ^3k  “  ^k’  ^4k  "  Pk’  ^5k  “  ^6k  " 


Solving  for  from  Eqs  22c  in  terms  of  one  finds 


r(^)  -  a 

S  “  \s^s 


(k,s  =  1,2, . . .6) . 


Substituting  Eq  25  into  Eq  22a  gives 


t  I  s  ,  ^/^jk  <^2^  ^  ®jk^  2^  ^(k+3)s  ^  2^^  ^ 

3—1  k— 1 


Equations  22b  and  26  constitute  a  system  of  homogeneous  linear  algebraic 
equations  in  The  existence  of  a  nontrival  solution  for  re- 

quires  the  vanishing  of  the  coefficient  determinant 

|a(6)I  =  0,  (27) 

where  A (6)  is  a  six  by  six  matrix  involving  6  in  a  transcendental  form.  Thus, 
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Eq  27  is  a  transcendental  characteristic  equation  for  the  standard  eigenvalue 
problem.  It  has  a  very  complicated  structure  as  can  be  seen  from  the  coeffi¬ 
cients  of  in  Eqs  22  and  26,  and  the  detailed  expression  for  A  (5)  is  not 

given  here.  The  investigation  of  the  characteristic  equation  requires  the 
employment  of  standard  numerical  techniques  such  as  the  Muller's  method  [28] 
with  the  aid  of  a  digital  computer.  The  eigenvalues  6^  obtained  from  the 
numerical  solution  of  Eq  27  give  important  information  concerning  the  behavior 
of  the  edge  stresses  and  displacements.  Due  to  the  positive  definiteness  of 
strain  energy  of  the  elastic  body  and  the  condition  in  Eq.  18,  the  eigen¬ 
value  of  6  bounded  by 
n 

-1  <  Re[6  ]  <  0  (28) 

n 

characterizes  the  order  of  the  inherent  singularity  of  the  boundary— layer  or 
free-edge  stresses  in  a  composite  laminate.  Thus,  for  small  values  of  r,  the 
asymptotic  stresses  are  proportional  to  ^  provided  that  5^  satisfies 

Eq  28. 
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4,  DEGENERATED  CASES  -  CROSS-PLY  COMPOSITE  LAMINATES 


In  the  case  of  a  cross-ply  composite  laminate,  i.e.,  laminae  with  0°  and 
90°  fiber  orientations  only,  the  in-plane  stresses  and  the  out-of-the-plane 
shear  stresses  are  uncoupled  by  virtue  of  the  material  symmetry  in  each  lamina. 
We  shall  concentrate  our  study  on  the  four  stress  components  of  interest  here, 
CT^,  general  expressions  for  displacements  and  may  be 

simplified  as 


u  =  -  +  UCx,y)  + 


(29a) 


V  =  -  j  +  V(x,y)  +  m^x  -  +  v^. 


(29b) 


w  =  (A^x  +  A^y  +  A2)S22Z  +  Wj^y  -  co^x  +  w^. 


(29c) 


“z  '  V-  '^2='  *  L2“y>'®33- 


(29d) 


Following  the  same  procedure  shown  in  the  previous  section,  the  govern¬ 
ing  partial  differential  equations  are  uncoupled  and  may  be  written  as 


L^F(x,y)  =  0, 


(30a) 


where  L,  is  defined  as  before 
4 


~  ^22  ^^^12  ^66^9x^8y^  ^11  37’’" 


(30b) 


The  homogeneous  solution  for  Eq  30a  may  be  obtained  in  a  simplier  form, 


F(x,y)  =  \  F  (x  +  P.  y). 


where  y,  are  the  roots  of  the  follow’ing  algebraic  equation: 


£^(y)  =  S^^y^^  +  C2S^2  ^22  " 
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The  homogeneous  stress  and  displacement  solutions  are  then  given  as 


(33) 

(34) 

(35) 


We  shall  choose  the  form  of  Fj^(x,y)  as 


F^(\)  =  I  C^zJ'^V[(5+l)(6+2)] 

k=l 


where  C,  and  6  are,  as  before,  arbitrary  complex  constants  to  be  determined 

tC 

later.  Imposing  the  homogeneous  traction-free  boundary  conditions  along  the 
free  edges  and  the  continuity  condition  along  the  ply  interface,  one  can 
precede  with  the  same  procedure  outlined  in  the  previous  section.  Then,  the 
eigenvalues  and  eigenfunctions  can  be  determined  in  a  manner  similar  to  those 
in  the  previous  cases. 


5.  NUMERICAL  EXAMPLES 


From  the  structure  of  the  governing  partial  differential  equations  and 
the  homogeneous  solution  for  the  problem,  it  is  clearly  seen  that  the  asymp¬ 
totic  stress  and  strain  fields  in  the  vicinity  of  the  edge  are  governed  by 
the  singular  terms  with  the  strength  of  stress  singularity  5^  determined 
from  the  eigenvalue  analysis.  Examining  the  structure  of  Eqs  22b  and  26,  it 
is  obvious  that  the  eigenvalue  solutions,  and  therefore,  the  edge  stress 
singularities,  are  related  to  laminar  constitutive  properties  and  fiber  ori¬ 
entations  of  adjacent  plies. 

Consider  a  composite  laminate  with  ply  properties  typical  of  those  used  in 

earlier  studies  [5]  (values  of  Gr£^=1.5G^'j<  and  .85  Vlt  have  been  found  to 

make  only  a  few  %  difference  in  6n) : 

E  =  20  X  10®  psi,  E„  =  E  =  2.1  X  psi, 

iu  1  Z 

Glt  =  =  0.85  X  io6  psi,  (37) 


0.21, 


where  the  subscripts,  L,  T,  and  z  refer  to  the  fiber,  transverse,  and  thick¬ 
ness  directions  of  an  individual  ply,  respectively.  The  influence  of  material 
properties  of  composite  plies  on  the  boundary  layer  stresses  may  be  related  to 
the  roots  of  the  characteristic  equation,  Eq  12a.  With  the  lamina  proper¬ 
ties  given  above,  the  roots  of  the  characteristic  equation  for  the  graphite- 
epoxy  lamina  of  different  fiber  orientations  0  are  shown  in  Table  1.  It 
appears  that  all  the  six  roots  y^^  are  purely  imaginary  by  virtue  of  the 
material  properties  of  Eqs  37.  Furthermore,  the  y^^  for  +0  ply  is  the 
same  as  those  for  due  to  the  in-plane  rotation  of  fiber  directions. 


Based  on  the  material  constants,  y^,  p^,  q^  and  t^  obtained  for  the 
graphite-epoxy-,  the  transcendental  characteristic  equation,  Eq  27,  can 
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Table  1 


Roots  of  Characteristic  Equations  for  Graphite-Epoxy 
Composite  System  with  Fiber  Orientation  0 


±0 

•^1,2 

’'3,4 

^5,6 

15° 

±0.88782 

i 

±1.10301 

i 

±1.56776 

30° 

±0.86222 

i 

±1.06956 

i 

±2.53630 

45° 

±0.80902 

i 

±1.03503 

i 

±3.44438 

60° 

±0.73316 

i 

±1.01323 

i 

±4.15870 

O 

in 

±0.66324 

i 

±1.00294 

i 

±4.61201 

*0is  the  angle  measured  counterclockwise  from  the  positive 
z-axis  to  the  fiber  direction 
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be  solved  numerically  to  provide  the  eigenvalues  for  the  homogeneous 
solution.  For  illustrative  purposes,  the  first  twelve  non- integer  eigen¬ 
values  associated  with  the  stress  solution  for  the  free  edge  of  a  45‘’/-45° 
graphite-epoxy  are  shown  in  Table  2.  Eigenvalues  6^  smaller  than  -1  are  ex¬ 
cluded  for  the  reasons  given  in  the  previous  section.  It  is  seen  that  there 
exists  one  and  only  one  eigenvalue  which  satisfies  the  required  constraint 
condition, Eq  28, for  this  case,  i.e.,  6^  =  -0.02557.  In  fact,  the  eigenvalue 
5^  is  the  strength  or  order  of  the  free-edge  or  boundary-layer  stress  singu¬ 
larity,  which  is  of  major  concern  in  this  study.  The  fact  that  there  is 
only  one  5^  which  meets  Eq  28  is  observed  in  all  cases  with  various  other  fiber 
orientations  studied;  the  only  difference  is  that,  for  each  case,  5^  possesses 
a  different  value.  Higher-order  eigenvalues  occurring  as  integers  (including 
zero)  and  as  complex  conjugates  always  exist  and  should  be  included  for  de¬ 
termining  the  complete  solution  when  remote  boundary  conditions  are  matched 
by  a  numerical  method  to  be  discussed  in  an  associated  report  [29]. 

For  the  commonly  used  (±0)  angle-ply  graphite-epoxy  composite,  as  anti¬ 
cipated,  the  order  of  the  boundary— layer  stress  singularity  is  a  function  of 
the  fiber  orientation  0.  Numerical  results  of  6^  for  each  of  the  ±9  fiber 
orientations  are  calculated  and  shown  in  a  graphic  form  in  Fig.  2.  It  is 
clearly  seen  from  the  Figure  that  the  composite  free  edge  associated  with  the 
laminate  of  an  approximately  (±51°)  ply  orientation  possesses  the  strongest 
boundary-layer  stress  singularity.  As  the  0  changes  to  either  directions, 
the  order  of  the  stress  singularity  6^  decreases  rapidly.  Its  value  con¬ 
verges  to  zero  for  the  cases  of  0  =  0°  or  90°,  since  the  two  plies  become 
identical  with  orthotropic  elastic  properties. 
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Table  2 

“fc 

First  Twelve  Non-Integer  Eigenvalues  for  Free-Edge  Stress 
Solutions  in  (±45®)  Graphite-Epoxy  Composite 

-2.5575658  E-2 

8.8147184  E-1  ±  i  2.3400497  E-1 
1.5115263  E  0  ±  i  7.9281732  E-1 

2.3389433  E  0  ±  i  1.1158402  E  0 

3.0913532  E  0  ±  i  1.7360464  E  0 

3.9520023  E  0  ±  i  2.0287146  E  0 

4.7440929  E  0  ±  i  2.5683871  E  0 

5.6021457  E  0  ±  i  2.8588510  E  0 

6.3962635  E  0  ±  i  3.3652707  E  0 

7.2565174  E  0  ±  i  3.6575937  E  0 

8.0497237  E  0  ±  i  4.1479983  E  0 

8.9120567  E  0  ±  i  4.4407609  E  0 

ie 

Integers,  0,l,2,...n,  are  always  eigenvalues  obtained  from  Eq  27 
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In  the  case  of  a  free  edge  associated  with  two  plies  of  more  general  fiber 
orientations,  instead  of  the  sjnnmetric  +Q/-0  configuration,  solutions  for  the 
eigenvalues  5^  are  obtained  also.  To  illustrate  the  nature  of  the  eigenvalues 
for  this  case,  we  examine  the  free  edge  associated  with  30°/®  fiber  orienta¬ 
tions  in  a  graphite-epoxy  composite,  where  ©  varies  from  7.5°  to  82.5°.  The 
first  few  non- integer  eigenvalues  for  various  0's  in  these  cases  are  given  in 
Table  3.  The  integers  (including  zero)  are  also  eigenvalues,  but  not  included 
in  the  Table.  The  case  of  30°/ 30°  graphite-epoxy  is  not  included  either  since 
the  two  plies  are  identical.  Again,  it  is  observed  from  the  Table  that  there 
exists  only  one  6^  which  meets  the  requirement  of  Eq  28  and  gives  the  dominant 
singular  stress  state  at  the  edge  of  each  of  the  30°/©  composite  laminate. 

The  degenerated  case  of  cross-ply  composite  laminates  which  are  discussed 
in  the  previous  section  is  investigated  also.  The  eigenvalues  for 
the  free-edge  stresses  in  a  graphite-epoxy  composite  with  0°/90°  lamination 
are  given  in  Table  4.  The  dominant  stress  singularity  in  the  present  case  has 
an  order  of  magnitude  similar  to  those  in  (±0)  angle-ply  composites  and  in  more 
general  (0^/02)  laminates.  It  is  noted  that  the  orders  of  the  boundary-layer 
stress  singularity  for  both  angle-ply  and  cross-ply  composites  are  generally 
much  weaker  than  those  associated  with  other  typical  elastostatic  singular 
stress  problems  such  as  elastic  crack  problems.  The  relatively  weak  singularity 
for  the  boundary-layer  stresses  Introduces  some  unique  features  as  well  as 
difficulties  for  the  evaluation  of  the  boundary- layer  effects  in  composite 
laminates,  which  will  be  discussed  in  an  associated  report  [29]. 
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First  Five  Non-Integer  Eigenvalues  for  Free-Edge  Stresses 
Associated  with  (SO®/©)  Graphite-Epoxy  Composite 
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Integers,  0,1,2,.. ..n,  are  also  eigenvalues. 


Table  4 


First  Twelve  Non-Integer  Eigenvalues  for  Free-Edge  Stresses 
in  Cross-Ply  Graphite-EpoKy  Composite* 


-3.33888  E-2 
8.80268  E-1 

1.41674  E  0  ±  i  3.93303  E-1 

1.65345  E  0  ±  i  6.85523  E-1 

2.83449  E  0  ±  i  1.76219  E  0 

3.75294  E  0  ±  i  1.1853E  E  0 

4.29235  E  0  ±  i  2.66884  E  0 

5.70726  E  0  ±  i  3.57190  E  0 

5.79010  E  0  ±  i  1.52461  E  0 

7.12293  E  0  ±  i  4.48145  E  0 

7.81068  E  0  ±  i  1.76401  E  0 

Integers ,  0,1,2,3...,  are  also  eigenvalues 
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6.  SUMMARY  AND  CONCLUSIONS 


A  study  of  boundary-layer  stress  singularity  in  both  angle-ply  and  cross- 
ply  composite  laminates  has  been  presented.  Formulation  of  the  problem  is 
based  on  Lekhnitskii’ s  complex-variable  stress  functions  and  basic  relation¬ 
ships  in  the  anisotropic  elasticity  theory.  An  eigenfunction  expansion  method 
has  been  developed  to  obtain  the  homogeneous  solution  for  the  coupled  govern¬ 
ing  partial  differential  equations  for  the  problem.  Angle-ply  and  cross- 
ply  composites  as  well  as  more  general  laminates  have  been  studied.  The 
strength  of  boundary-layer  stress  singularity  for  each  case  has  been  determined 
to  illustrate  the  fundamental  nature  of  the  edge  effects  in  composite  materials . 

Based  on  the  information  obtained,  the  following  conclusions  may  be  drawn: 

1.  Boundary-layer  or  free-edge  stress  field  in  a  composite  laminate 
is  inherently  singular  in  nature  due  to  the  geometric  and  material 
discontinuities . 

2.  The  order  of  boundary-layer  stress  singularity  can  be  determined 
by  solving  for  the  transcendental  characteristic  equation  ob¬ 
tained  from  the  homogeneous  solution  of  the  governing  partial 
differential  equations . 

3.  The  boundary-layer  stress  singularity  depends  only  upon  material's 
elastic  constants  and  fiber  orientations  of  adjacent  plies  in 
composite  laminates. 

4.  For  angle-ply  and  cross-ply  composites  as  well  as  more  general 
laminates  the  order  of  boundary-layer  stress  singularity  is  very 
weak  in  general.  In  a  graphite-epoxy  system,  for  example,  5^  is 
much  smaller  than  other  kind  of  singular  stress  problems  in 
elastostatics  such  as  elastic  crack  problems. 
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APPENDIX  1 


Expressions  for  H^^(0) 


in  Equations  19  and  20 


H, 


Ik 


H 


2k 


H, 


3k 


H 


4k 


H 


5k 


6k 


H 


7k 


H 


8k 


(Pj^sinB  +  cos6)^ 

-rikCl-'k®^’^®  +  COS0) 

-(lJ^sin0  +  cos6)  (Pj^cos0  -  sin0) 
(yj^cos0  -  sin6)^ 
rik(MkCOS0  -  sin0) 
p^cos0  +  qj^sin0 

-Pj^sin6  +  qj^cos6 
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